Critical Domain Walls in Ashkin- Teller Model 
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We introduce a new family of critical curves for the Ashkin- Teller model on the critical line. 
Starting from known results for some special points, the fractal properties of these curves are in- 
vestigated. A relation cos 2 (27rD) = 4 1 2 ^ s proposed, valid all along the critical line, relating the 
thermal exponent v and the fractal dimension of cluster boundaries D. High-precision Monte Carlo 
results on interfaces generated in the triangular lattice are presented in support of the conjecture. 
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INTRODUCTION 

Studying critical systems has been in the interest of 
physicists for at least five decades. In two dimensions 
many exact methods to study critical systems were in- 
vented by physicists and mathematicians, such as the 
Yang-Baxter equation [3, conformal field theories (CFT) 

and recently Schramm-Loewner evolution (SLE) 
Most of the critical systems in two dimensions can be 
formulated in terms of fluctuating loops: for example, in 
the most familiar case, i. e. the Ising model, they simply 
correspond to domain walls between regions of opposite 
magnetization. The SLE is based on direct investigation 
of these critical loops and is based on probability tech- 
niques (for a review, see [H). According to SLE all of 
the conformally invariant curves in two dimensions can 
be parametrised by the SLE drift k, related to the fractal 
dimension of the curves via D = 1 + ^ . 

One of the most interesting statistical systems to study 
from the SLE point of view is the Ashkin- Teller (AT) 
model [5|. This model is interesting at least from two 
points of view: first, it has a rich phase diagram including 
a critical line, and, secondly, it has an interesting physical 
realisation as Selenium adsorbed on the Ni(100) surface 
Q. Following the above motivations we investigate the 
critical loops in this model on the critical line. 

The definition of the Ashkin- Teller model on an ar- 
bitrary graph is as follows: on each vertex r of the 
graph lives a field s r with values s r — e l ^ q , where 
q = 0, 1, 2, 3. Then the partition function of the model is 
z = E{ Sr } Ilr.r' W(s r , Sr'), where: 



W(s r , S r t) = (l + XlS*S r ' + X2st 2 s 2 



I, +x 1 s* 3 sl,) , (1) 



with the product over pairs of adjacent sites. The above 
partition function is apparently Z4 symmetric and for 
x\ = x-i reduces to the four-state Potts model. The 
model was solved exactly on the square lattice by map- 
ping it to the 6-vertex model and some of the critical 
exponents were found bymapping the model to Solid-on- 
Solid (SOS) model (0,1,11; see also \\J$ for the phase 
diagram of the anisotropic case). The self-dual Ashkin- 
Teller model on the square lattice is described by the line 



X2+2xi = 1; moreover, it is exactly solvable all along this 
line. Some special points along this line are in the univer- 
sality class of well known models: the point x\ = X2 = 3 
is in the universality class of the four-state Potts model. 
The model is critical for all points on the self-dual line 



with Xi > i. At the point X\ 



sin(ir j 16) 



called Fateev- 



3 . sm(37r/16) : 

Zamolodchikov (FZ) point, the model is fully integrable 



and can be described by Z4 parafcrmionic CFT llj. By 
an easy change of variables, see Eq.[?l one can show that 
the above partition function corresponds to the Hamil- 
tonian of two coupled Ising models which decouple for 
X2 = x\, so for x\ = 1/2 — 1 we have the Ising universal- 
ity class. Finally, at xi = 1+2 g m g the system falls into 
the universality class of the XY model. Since the critical 
properties of the fluctuating curves at the above special 
points on the critical line are well-known, one can use 
them to extrapolate the results to other points. 

From the CFT point of view it is well-known that the 
field theory describing the Ashkin- Teller model at the 
critical line is the c = 1 orbifold conformal field theory 
[HI, EH, EH- This CFT comes from compactifying the 
free bosonic field theory to a circle with radius r and 
then requiring a Z2 symmetry for the bosonic fields (for 



3, 13). to 



a review of this conformal field theory see 
fix the notation take the free field theory with action 
S = j- J d<f>d(j> and then compactify the bosonic field 
on a circle of radius r. The action is invariant under 
4> — * —<p. To get an orbifold CFT we should project all 
the operator content of the circle theory to the operators 
which respect this symmetry, however, this theory cannot 
be modular invariant without introducing some twisted 
operators which come from the discrete Z2 symmetry of 
the model. In the case of Z2 orbifold CFT we have four 
operators, two of them have conformal weights jg, and 
we will call them cr, which is reminiscent of the spin op- 
erator in Ising model. The other two have conformal 
weights jq and are called t. This theory is conformally 
invariant for all real values of r and by changing the ra- 
dius we can change continuously the critical exponents of 
the model except the twist operators which are invariant 
under a change of orbifold radius. On the square lattice 
we have the following equality between the radius and 
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the coupling of Ashkin-Thcllcr model on the critical line 
MM: 
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(2) 



This equation is valid only on the square lattice: one 
lattice-invariant way to write the above equation, which 
is also useful in numerical calculations, is: 



i-2-» 

v r A 



(3) 



where v is the thermal exponent. Some well-known 
points are the following: r — 2 is the four-state Potts 
model with 11 primary fields as a modular invariant orb- 
ifold conformal field theory [13]. r = a/3 describes Z4 
parafermionic CFT with 10 primary operators, which 
give s again a modular invariant partition function (see 
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17(1). r = V2 is in the Ising universality class, indeed 
it is a pair of decoupled Ising models. The interesting 
point about all these models is that we can write the 
compatification radius as r = ^/p, where p is an integer 
number: the CFTs corresponding to these orbifold field 
theories (studied in detail in [17|) have p + 7 primary 
operators. The equation r = ^Jp is reminiscent of the 
loop weights n in the Fortuin-Kasteleyn representation 
of the Q-state Potts model, i. e. n = \/Q. The compacti- 
fication radius on these particular points is also equal to 
the quantum dimension of the spin operator a on these 
points [H. 

There are different possibilities to define critical curves 
in the Ashkin- Teller model. In the spin representation 
one can think about the domain walls between one spin 
and the other three (see [2, or the domain walls 

between two definite spins and the other two. It is ob- 
vious that at the Ising point the latter ones, if properly 
chosen, give the domain walls of the Ising model, so we 
will focus on a "symmetric" choice recovering standard 
Ising interfaces on one of the two underlying Z2 systems. 
From well-known results, the SLE drift at the Ising point 
is k = 3. At the FZ point the critical curves are in the 
universality class of the 3-state Potts model 2l| and the 



SLE drift is k — -y . Finally at the four-state Potts model 



the SLE drift is k = 4. 
in a compact form as 



One can write the above results 
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If the fractal properties of loops on the critical line change 
continuously then the above simple function is the most 
promising function to describe the critical curves at the 
generic point on the critical line. The main claim of this 
paper is that Eq. (j4|) is valid all along the critical line. We 
numerically check this equality on the region between the 
Ising point and the four-state Potts model point. From 
the connection between boundary CFT and SLE, it is a 



known result that for all the three benchmark points the 
corresponding boundary changing operator is the energy 
operator e with conformal weight h e = We believe 
that the same operator is responsible for changing BC 
at the generic point on the critical line. To check the 
validity of the conjecture, we calculate numerically the 
fractal dimension of the curves and the thermal exponent 
v along the critical line, which are expected to satisfy: 



cos 2 (2ttL>) = cos 2 27T (1 
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DEFINITIONS, PROCEDURE AND 
ALGORITHM SETTINGS 



The numerical part of this work is conveniently ex- 
pressed in the coupling space /3, a, where: 

z = E ; («) 

{<T r ;-r r } r,r' 

S rr i = /3((7 r 0Y< + T r r r i ) + a(a r a r iT r T r i) . (7) 

Here, a and r are two Ising variables; the correspondence 
is s r = e ^ (py +ir r ), so that: 
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e 2f3 + e -2/3 + 2e -2a 



X2 



=2/3 



-2/3 



2e~ 2Q ' 
(8) 

In this language the Ising model has a = 0, and the 
Potts line is a = (3: the update algorithm used is an 
adapted version of the Swendscn-Wang prescription, al- 
ternatively applied to the two sublattices 22j. The nu- 
merical work was carried on in the triangular lattice be- 
cause in this case an interface can be defined in a unique 
and natural way (the dual lattice possesses only three- 
links joints). To validate Eq. [5] with Monte Carlo data, 
two steps are necessary: first, identify some points along 
the triangular-lattice critical line and numerically extract 
the associated exponent v\ then, measure the interface 
length S(L) at those points for a variety of system sizes 
and look for the fractal dimension as the exponent in the 
relation, valid at large L, 



S(L) = aL D 



(9) 



To locate a point and the associated v on the critical 
line, we measured, as in [23j | . the (density of) susceptivity 
for the operator S rr >, i. e. (x)/3;L = ((S - (S) ;L ) 2 ) 0;L ■ 
Measurements were taken for various periodic-boundaries 
finite sizes L, along a segment in the phase space sup- 
posed to cross the critical line. Since the power-law scal- 
ing form of this quantity near criticality is known, is it 
possible to numerically tune (3 C and v to the best collaps- 
ing of the curves onto the universal rescaled function: 



Fl\ x — (P ~ flc)L») = L 4 ~»(x)/3 ; i 



(10) 
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P 


a 






D(i/= xp ) 


A 


0.19356(8) 


0.1343 (3) 


0.747(5) 


1.418(1) 


1.418(2) 


B 


0.21694(1) 


0.09246(9) 


0.807(5) 


1.400(1) 


1.400(2) 


C 


0.24210(2) 


0.05051(8) 


0.858(4) 


1.391(1) 


1.391(1) 



TABLE I: The couplings, thermal exponents and fractal di- 
mensions for the three points A, B and C. ^ exp is the numerical 
estimate from the data-collapsing fit. The fractal dimensions 
are the measured one and that coming from u cxp via Eq. [5] 
Statistical errors (in parenthesis) affect the last digits. 
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Universal rescaled curves for A,B,C 
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S(L), points A,B,C 



For the direct extraction of the fractal dimension at 
a given critical j3, a point, an interface must be cre- 
ated, via fixed boundaries, to run across the system (as- 
sumed of rectangular shape). This is achieved by forcing 
cr = +( — ) on the left (right) half of the system's bor- 
der, neglecting the value of r. This implements, on the 
tr layer, a vertical interface of the kind defined above. 
Given a configuration, the associated interface length S 
is obtained by rccolouring isolated sign-clusters to leave 
only two regions, and then counting the links connecting 
spins of opposite sign [lj)[ . Note that sometimes part of 
the interface can lie on the border (for a total of Sb(L) 
steps): this phenomenon is expected to be negligible for 
large systems. 
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FIG. 1: Top: the rescaled susceptibility curves F(x) for points 
A, B, and C (bottom to top). Bottom: log- log plot of the data 
for S(L) at the A, B, C points. The lines are the best fits to 
EqM 



RESULTS 

On the triangular-lattice critical line, three points la- 
beled A, B and C, evenly spaced between the Potts 
(af = /3 C P = *2tp) and the Ising = ^) points, were 
chosen to test the conjecture. The critical couplings and 
v indices were extracted as described above, by sampling 
the susceptibility over a range of ~ 30 couplings at 6 
square-shaped system sizes up to L — 400 and with a 
statistics of 80.000 configurations. The data-collapsing 
procedure worked well, allowing for a precise identifica- 
tion of the A, B, C points and their index v (Fig. [TJ top). 

At the three numerically found points, interface 
lengths S(L) were measured over 325.000 configurations, 
for 40 different square-shaped systems of size L from 120 
to 900 (Fig. [TJ bottom). To compare the collected data 
with Eq. [9j we progressively discarded small sizes up to 
finding an acceptable \ 2 ■ This typically happened at 
about L min ~ 400, and provided good plateaux for es- 
timating D with some accuracy (Fig. [5J . The obtained 
Monte Carlo results, Table [I] can be used to test Eq. [5] 
Fig. [3] clearly supports our conjecture. 

To guarantee the results are not biased by the verti- 
cal walls of the system, a separate test was conducted at 
the Ising point on systems with horizontal extent 1, 2, 4, 8 
times the vertical one: all the D extracted were compat- 
ible, provided the measurements arc taken with L large 



Plateaux for D, points A,B,C 



1.39 
1.38 



point A i — i 1 

point B « 

Point C • r-f~H- ~ }-}-- [- 












if i 11 

****f ! 







150 200 250 300 350 400 450 500 550 



L_min 

FIG. 2: The fitted values of D approach a plateau at large 
enough L m i n . Lines represent D exp . 



enough. 

Finite- L systematics seem to be of concern especially 
near the Ising point: indeed, tests performed at the 
FZ and Ising points (the former being very close to A), 
while proving the overall consistency of the method used, 
showed that the latter suffers from some finite-size bias 
(of about 0.4% on the final D). Nevertheless, the rela- 
tion between the measured v and D is preserved: the 
finiteness of the systems used induces a sort of "effective 
thermal exponent" , which enters both measurements in 
the same way. 

It has been verified, furthermore, that the contribu- 
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Fractal dimension, points A,B,C 




fruitful discussions. M. A. R. thanks M. Ghascminczhad- 
haghighi for early discussion of this problem. 



FIG. 3: Comparison between the measured fractal dimension 
D in A, B, C (points) and the estimate coming from Eq. [5] 
(outer lines mark the one-sigma compatibility strip). 



tion of the bordering part of the interface effectively 
vanishes for large systems: numerical results support 
the paramctrisation Sb{L) = s + bL n , with a small 
(~ 50) additive constant and exponent never exceeding 
half the fractal dimension: the ratio — has the values 

n 

2.95(1), 2.46(1), 2.11(1) for A, B, and C respectively (and 
is compatible with 2 at the Ising point). 



CONCLUSIONS 

In this Paper we proposed a generalisation of the rela- 
tion between the SLE drift of critical spin-cluster inter- 
faces in the AT model and the compactification radius 
of the corresponding orbifold CFT: starting from known 
results valid for some special points, we claim the result 
to be valid all along the critical manifold. High-precision 
Monte Carlo measurements support this conjecture. The 
validation is enabled by the possibility of relating the ra- 
dius r to the thermal exponent v of the statistical model: 
however, the measurements of v and D come from two 
conceptually different frameworks, so that the agreement 
between the two quantities is highly non-trivial. 

The Monte Carlo fractal dimension measurements were 
conducted on those points resulting effectively critical at 
the system sizes considered: although they do not lie ex- 
actly on the infinite- volume critical line [24j . it would be 
incorrect to compare the v measured there and the D 
extracted from simulations displaced in the phase space 
(especially because the exact physical flow to criticality 
is unknown and this would introduce errors). Further- 
more, it would be desirable to check the stability of these 
results within a different discretisation scheme: measure- 
ments were indeed taken also on the square lattice, gen- 
erally confirming Eq. [5j but there the very definition of 
an interface has some ambiguities pEj , associated with a 
"four-lines vertex" . A detailed investigation on the fea- 
tures of the square lattice will appear elsewhere (2(| . 
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